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If 7-ray bursts (GRBs) are accompanied by gravitational wave bursts (GWBs) the correlated output 
of two gravitational wave detectors evaluated in the moments just prior to a GRB will differ from 
that evaluated at times not associated with a GRB. We can test for this difference independently of 
any model of the GWB signal waveform. If we invoke a model for the GRB source population and 
G WB radiation spectral density we can find a confidence interval or upper limit on the root-mean- 
square GWB signal amplitude in the detector waveband. To illustrate we adopt a simple, physically 
motivated model and estimate that initial LIGO detector observations coincident with 1000 GRBs 
could lead us to exclude, with 95% confidence, associated GWBs with Tihms ^ 1-7 X 10~ 22 . This 
result does not require the detector noise be Gaussian or that any inter-detector correlated noise be 
measured or measurable; it does not require advanced or a priori knowledge of the source waveform; 
and the limits obtained on the wave-strength improve with the number of observed GRBs. 
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Gamma Ray Bursts (GRBs), which are known to lie 
at cosmological distances, likely arise from shocks in a 
relativistic fireball that is triggered by rapid accretion 
on to a newly formed black hole The violent for- 
mation of a black hole is likely to produce a substantial 
gravitational wave burst (GWB); thus, we expect GRBs 
to be preceded by GWBs that are candidate sources for 
the new gravitational wave (GW) detectors now under 
construction ||||. 

Proposed GRB progenitors include coalescing or merg- 
ing binary systems {e.g., neutron star or black hole bina- 
ries, white dwarf - black hole mergers, Helium core - black 
hole mergers) and hypernovae or collapsars {e.g., failed 
SNe lb and single or binary Wolf-Rayet star collapse) . 
In fact, statistical evidence points to at least three dif- 
ferent subclasses of GRBs [| ; so, the actual progenitors 
may include these as well as other systems. Matched 
filtering (MF) — the method at the focus of most of 
the gravitational wave detection literature — requires 
detailed knowledge of the actual GWB waveform; con- 
sequently, it cannot be used to detect the distinct GWB 
associated with a GRB. Additionally, since GRBs occur 
at cosmological distances {i.e., z > 1), the signal-to-noise 
ratio (SNR) of any individual GWB burst will almost cer- 
tainly be insufficient for a high confidence detection with 
these new detectors. Detection techniques other than 
MF that aim to detect distinct GWBs will perform even 
worse. 

Here we suggest an alternative method for detecting 
a GWB / GRB association. If GWBs are associated with 
GRBs, the correlated output of two GW detectors will be 
different in the moments immediately preceding a GRB 
{on-source) than at other times, not associated with a 
GRB {off-source). A statistically significant difference 



between on- and off-source cross-correlations would sup- 
port a GWB/GRB association and, by implication, rep- 
resent a detection of gravitational waves by the detector 
pair. (While we focus on GRBs in this paper any plausi- 
ble class of astronomical events can serve as a trigger.) 

We can test for a difference between the on- and off- 
source cross-correlation populations using Student's t- 
test, without specifying an a priori model for the signal 
waveform, source, or source population. The observed t 
statistic can also be used to set a confidence interval (CI) 
or upper limit (UL) on the root-mean-square (RMS) am- 
plitude of GWB signals associated with GRBs, where the 
average is over the source population. If we specify even 
a rudimentary GWB source model, the CI/UL constrains 
the model. 

We restrict attention here to the two full-length LIGO 
detectors (denoted 2>;, i = 1,2). These detectors are 
nearly identically oriented and lie ~ 3000 Km apart. We 
do not require that the detector noise statistics be Gaus- 
sian, although we do require here that they be approxi- 
mately stationary. Without loss of generality we assume 
the noise has zero mean and denote its one-sided power 
spectral density (PSD) by Si(f). Noise cross-correlated 
between T>\ and T>i is assumed only to be stationary (but 
not necessarily Gaussian) and weak compared to the in- 
trinsic noise of both the detectors. 

a. The on-source and off-source distributions. Sup- 
pose that a GWB, associated with a GRB, is incident 



from direction ft on the GW detector T>i at time t a 



(i) 



The lag 5t, equal to t { a' — 1^\ depends only on n, which 
we know from the GRB observation. The lag is also the 



same as the difference t 



(2) Al) 



t\ x> , where t-fi is the arrival 



time at detector T>i of the GRB. 

Assuming that any associated GWB precedes the 
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GRB, focus attention on the output Xi(t), for < 
*7 — t < T, of detector T>i. Choose T as long, but 
no longer, than necessary to insure that Xi includes the 
possible GWB signal. From the 2^(4) we compute the 
weighted cross-correlation 



X := (xi,x 2 ) 



dtdt'xx{t^ -t)Q(\t-t'\)x 2 (t^ -t'). (1) 



The filter kernel Q is at our disposal: we discuss its choice 
in section [] below. 

The collection of X computed for each of N on GRBs 
form a set X on of on-source events. To complement that 
set, we also construct a set X g of N a g off-source events, 
using data segments Xi corresponding to random sky di- 
rections and arrival times not associated with any GRB. 

The sample sets X a g and X on are drawn from popu- 
lations whose distributions we denote p Q s and p on . For 
T much greater than the detector noise auto- and cross- 
correlation times, the central limit theorem implies 
that p ff is normal, 



p oS (X) ~N(X; iioB,a), 

with mean and variance 

H oS := E[(ni,n 2 )], 
a 2 := E ((ni,ra 2 ) - ^ Q s) 2 



(2a) 

(2b) 
(2c) 



Here rii(t) denotes noise from detector i and E[-] rep- 
resents an ensemble average across the detector output. 
Note that /x g- is just the detector noise cross-correlation 
evaluated at the lag St. (A weak assumption behind equa- 
tion ^a| is that the noise cross correlation does not vary 
significantly over the lag St, which should be the case for 
terrestrial noise sources.) 

Now suppose that GRBs are preceded by GWBs. El- 
ements of X on then take the form 

X = (rii,n 2 ) + (hi,n 2 ) + (ni,h 2 ) + (hi, ha) , (3) 

where hi(t) is detector i's response to the incident GWB. 
Define Pi by 



Pi- 4/ df\hi(f)\ 2 /Si(f). 
Jo 



(4) 



If Pi, the average of Pi over the source population, is 
much less than unity, then p on is also a normal distribu- 
tion with variance a 2 and mean 



^on — MofF + S, 

s := (hi,h 2 ) 



where 



(5a) 
(5b) 



and s is, again, an average of s over the source popula- 
tion. 



b. Student's t-test. Pose the null hypothesis 

Ho : Poff(X) = Po n(X) . 



(6) 



Rejecting H supports a GWB/GRB association. Since 
p on and p Q ff are normal and differ, if at all, only in their 
means, we can test H using Student's t-test Q. 
The t statistic is defined from X on and X a g by 



Aon - Moff / N on N ofi 



N 



N oS ' 



a 2 on (N on - 1) + a 2 oS (N oe - 1) 



N, 



off 



(7a) 
(7b) 



where p, on and p, Q g (a 2 n and a 2 s ) are the sample means 
(variances) of X on and X g, respectively. 

The expectation value of t, averaged over the source 
population and across the detector noise processes, is 



M, :=E[t] = lJ N °» N ° S . 
m U a V iVon + Non 



(8) 



The relative orientation of the two LIGO detectors guar- 
antees that hi(t) and h 2 (t) are very nearly identical. It 
follows that, for LIGO, s is non-negative; correspond- 
ingly, the expectation value E[t] of t is positive in pres- 
ence of a GWB / GRB association and zero otherwise. 

The actual value of t, given any observed sets X on and 
X g, will vary from \x t . The distribution of t is normal 
for large N OD + N a g and is, more generally, tabulated in 
any standard statistics text Consequently, we can 
find a to such that, when Hq is true (/i t = 0) t is greater 
than to in less than a fraction a (e.g., 5%) of all obser- 
vations. This is our test: if we observe t greater than to 
we reject Ho and conclude that we have found evidence 
of a GWB/GRB association with significance 1 — a (e.g., 
95%). 

c. The filter kernel Q. If we knew the signal hi(t) 
corresponding to each GRB trigger we could construct a 
Q that maximizes s: 



Q(t)= / dfe 



,2-rrifT 



hi(f)h* 2 (f) 

si(\f\)s 2 (\f\y 



(9) 



where hi is the Fourier transform of hi. For the LIGO 
detectors, the hi are identical. Denoting their common 
functional form h(t), the optimal Q depends only on the 
GWB signal spectral density \h(f)\ 2 . 

More generally, we can put any knowledge we have of 
the signal's spectral density shape into Q. For LIGO we 
can choose Q to match the signal model irrespective of 
the details of the possible GWB waveforms if \h(f)\ 2 is 
independent of signal parameters. This happens, for in- 
stance, in the case of the quadrupole formula waveform of 
an inspiraling binary. For GWBs associated with GRBs 
there is no reason to believe that \h(f)\ 2 will be known a 
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priori, let alone that it have this special property. Lack- 
ing detailed knowledge, we recommend adopting Q given 
by equation || with \h(f)\ 2 assumed to be unity in the 
detector bandwidth. 

d. Setting upper limits. Having specified Q we can 
test Hq as described above to rule on the presence of a 
GWB/GRB association, independently of any model for 
the GWB signal or its source. Alternatively, we can use 
the observed t to establish a confidence interval (CI) or 
upper limit (UL) on fj, t , and hence s, which is related to 
the GWB wave-strength (cf. eq. |, |t]). If we invoke a 
model for the spectral density |/i(/)| 2 and spatial distri- 
bution of GWB/GRB sources, this becomes a physical 
constraint on the model. 

To obtain a CI/UL on s given an observed t, we con- 
struct a confidence belt 0, of desired confidence, in the 
t-fi t plane. We follow the construction of |Q, which uni- 
fies the treatment of CIs and ULs. ( (7[ table X] tab- 
ulate UL/CIs appropriate to our case: the mean /i of 
a Gaussian variable x when /i > and x itself is the 
observation.) From observed sets X on and X a g the corre- 
sponding t is computed. From t and the confidence belts 
the corresponding CI/UL on s can be read-off. 

To interpret the CI/UL, imagine an ensemble of obser- 
vations and corresponding pairs of sample sets X on and 
X ff. Corresponding to this ensemble of observations is 
an ensemble of values t and CI/ULs. If the confidence 
level is e (e.g., 95%), then a fraction e of these CI/ULs 
will include the actual value of s (and a fraction 1 — e will 
not). It is in this sense that the CI/UL corresponding to 
the observed t is said to bound s with confidence e. 

To measure the effectiveness of the proposed test we 
evaluate the UL most likely to be placed on s if Hq is, in 
fact, true. When Hq is true the most likely observed t is 
zero. Denoting the corresponding UL on /x t as /it, m ax the 
UL on s is 



• N, 



off 



N on N oS 



Mi, max y/2/Wy (AU = N oS = 
Mt, max /VAU (AW>AU) 



AL) 



(10a) 
(10b) 



Since the duty cycle of GRBs is low, the size of the off- 
source sample can be made much larger than the size of 
the on-source sample. Even if both sample sets are the 
same size, however, the limit obtained will be weaker by 
only a factor of 2 1 / 2 . 

The upper limit /it, m ax corresponding to an observed 
t of zero and different degrees of confidence is given in 
0, table X]. For reference we note that /it, max is 1-00 for 
68.27%, 1.64 for 90%, 1.96 for 95%, and 2.58 for 99% 
confidence. 

A derived CI/UL on s implies, within the context of a 
GWB/GRB source model, a CI/UL on the RMS GWB 
signal amplitude in the detector band, with the average 
over the source population. As an example, suppose that 



each GRB is accompanied by the formation of a several 
solar mass black hole and a corresponding millisecond 
timescale GWB in the source rest frame. Assume further 
that |/i(/)| 2 is approximately constant in the correspond- 
ing KHz bandwidth B s . (This is consistent with numer- 
ical models of supernova core collapse [|]|| and with the 
formation or ring-down of all but the most rapidly ro- 
tating solar mass black holes p0|.) At the detector, the 
signal power from a source at redshift z lies in the band- 
width B s /z' , where z' is equal to 1 + z. 

For simplicity, assume that the detector noise PSDs 
Si(f) are identical and equal to a constant So in the de- 
tector bandwidth Bd, which we take to be approximately 
100 Hz about a central frequency of 150 Hz. Outside the 
detector band we set Sj equal to infinity. (This is a rough 
approximation to the actual shape of the noise PSD of 
LIGO |11|.) Finally, note that B s is much larger than 
Bd, so that B s /z' completely overlaps Bd for some large 
range of z'. 

With these assumptions, 



df\h(f)\ 2 Q(f) = 



2A 2 B d 

°0 



T 
1 



#Si(|/|)S 2 (|/|)|Q(/)| 2 = 



and 
TB d 

2sr 



where A is defined by 



df Hf)\' 



2A 2 B, 



From equations O, [12 and O it follows that 



= E 



2y/2A 2 B d 



y/TBdSo 



2V2 A 2 B d 



y/TBdSo 



(11) 
(12) 

(13) 

(14) 



where we have replaced A 2 by its mean over the source 
population (a good approximation when A is sharply 
peaked about its mean.) From equations [lO] and [l4| and 
assuming that Hq is true we find 



A 2 < A 2 



Ih, 



2^2 



TB fl 



N n 



1/2 



~B~ d 



(15) 



with 7V ff 3> iV on and /it, max obtained from the confidence 
belt construction Q with t = 0. 

We expect that different GWBs will have different 
waveforms and durations. Define the RMS signal power 
in the detector band by 



h 



R.MS 



T Jf£B d 



dfWfW 



(16) 



where h(t) is the GWB waveform, r its duration in the 
detector band, and the average is over the source popu- 
lation. In our example — broadband bursts whose band- 
width includes the detector band — we can approximate 
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1/r by the detector bandwidth Bd- Combining equations 
|l6| , |l5| and [TJ] we find the UL on /irms : 



,2 < r 17x ln -221^ ft,max ( T 1000 \ 1/2 

^rms< [1-7x10 j 19g l 05giVo J 



So 



(3 X 10- 23 Hz-V2 



100 Hz 



3/2 



(17) 



The reference values of Bd and So are characteristic of the 
initial LIGO detectors @- For T (cf. eq. 0) we assume 
GRBs are generated by internal shocks in the fireball; 
then, the GRB/GWB delay is approximately 0.1 sec in 

To accommodate GRBs at 
0.5 sec. Finally, fit. max equal 



the source rest frame 12 
redshifts z < 4 we take T ^ 

to 1.96 corresponds to a 95% confidence UL 0. 

If, on the other hand, GRBs are generated when the 
fireball is incident on an external medium, then |l^, eq. 
3.6] with n t = 1, a = 1, E 51 = 10, and T > 100 gives 
a source rest-frame delay < 100 sec, in which case T 
should be 500 s and the corresponding UL on /irms is 
9.4 x 10" 22 . 

Two final notes are in order. To calculate the X (cf. 
eq. 0), which are at the heart of our analysis, we must 
know accurately the GRB source direction. Bright bursts 
in the BATSE3B catalog have positional accuracies of 
59 < 1.5° [jl4j]. The corresponding uncertainty in s is 
< 5%, which does not affect significantly the UL on s. 

Finally, the proposed BATSE follow-on — SWIFT - 
is not an all-sky GRB detector. It will have greater sen- 
sitivity than BATSE, but observe only a fraction of the 
sky at any one time. If SWIFT pointing favors the sky 
normal to the LIGO detector plane, LIGO's sensitivity 
to GWBs from observed GRBs will be maximized, in- 
creasing the sensitivity of the test described here. 

e. Conclusions. Gamma-ray bursts (GRBs), which 
are believed to be associated with the violent formation 
of a stellar mass black hole, may well be immediately pre- 
ceded by a gravitational wave burst (GWB). If we com- 
pare the correlated output of two gravitational wave de- 
tectors immediately preceding a GRB to the correlation 
at other times, not associated with a GRB, then a statis- 
tically significant difference is evidence for a GRB/GWB 
association. 

We can test for this difference — independent of any 
model of the GRB/GWB source or GWB waveform — 
using Student's t-test. Alternatively, we can set an up- 
per limit (UL) or confidence interval (CI) on the RMS 
GWB amplitude in the detector waveband, averaged 
over the source population. This CI/UL constrains any 
GRB/GWB model we do invoke. 

This analysis has several important advantages over 
matched filtering, which is the method at the focus of 
most of the gravitational wave detection literature. In 
particular, it becomes more sensitive as the number of 
observed GRBs increases, does not require any knowledge 



of the GWB waveforms, is insensitive to the presence of 
non-Gaussian detector noise, and does not require statis- 
tical independence of the detectors or knowledge of their 
correlated noise. It is thus a powerful addition to the 
growing arsenal of analysis techniques aimed at making 
gravitational wave detection an astronomical tool. 
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